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4.3 The Renormalization Group Methods 

Some typical results of MC simulations of a (i = 2 Ising ferromagnet are shown 
in Fig. 4.9. At high T (T — 2Tc), there is only short-range order, the spins 
form small clusters. The correlation length (approximately equal to the Unear size 
of the largest cluster) is small. Close (but above) Tc, somewhat larger patches 
in which most of the spins are lined up in the same direction begin to develop. 
When the system reaches Tc, these patches expand to infinite size, but fluctuations 
of smaller scale persist! At the critical temperature, spins form clusters at all 
lengthscales, including one infinite-size cluster (which cannot be seen on a finite 
system, but we know that the correlation length has to diverge at Tc). As a result, 
all scales of length must be included in a theoretical description. Another very 
important conclusion is reached if we substitute a block of 3 x 3 spins by a single 
spin with the orientation prescribed by the majority rule (block spin "up" if the 
majority of the spins in the block is "up", and vice versa). In this way the length 
scale of the lattice is changed by a factor 6 = 3 each time. This is a "real-space 
block-spin renormalization-group transformation". 

Examples of such transformations are shown in Figs. 4.10-4.12. Below Tc, 
the system has long-range order but, because the system is close to Tc, the fluc- 
tuations form large clusters of oppositely ordered spins. However, under each 
transformation the size of the unit cell increases, the fluctuations become smaller 
and smaller and those fluctuations, which are smaller in size than the lattice spac- 
ing, disappear. Under successive transformations, thus, the system becomes more 
and more ordered, it effectively approaches zero temperature, see Fig. 4.13. 

Above Tc (Fig. 4.11), there is no long-range order, spins form random "up" 
and "down" clusters. Under each transformation, the correlation length (expressed 
in units of the new unit cell) decreases, the clusters become smaller and smaller 
as if the temperature were higher, but they never disappear. At T = Tc (Fig 
4.12), the clusters are of all length scales, including one infinitely large cluster. 
The correlation length is infinite (exceeds the size of the lattice used in the MC 
simulation). Under successive transformations, it remains "infinite", there is no 
change in the distribution of cluster sizes, the system is invariant under scale trans- 
formations ! The system remains at the critical temperature and the critical point 
is di fixed point. Other (trivial) fixed points are T = and T = oo because the 
system away from Tc tends to one of these points under successive scale transfor- 
mations. Decrease of the correlation length under successive transformations is an 
extremely useful property because it makes the fluctuations uncorrected and we 
can solve the system using an approximate (e.g., mean-field) theory to calculate 
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Figure 4.9: Results of MC simulations of a ri = 2 Ising ferromagnet at different T. 
At high T (T = 2Tc; top), the correlation length is small. Close (but above) Tc 
(middle), somewhat larger clusters begin to develop. When the system reaches Tc 
(bottom), these clusters expand to infinite size, but fluctuations on smaller scale 
persist! (From: K.G. Wilson, Sci. Am. 241, 140 (1979).)) 
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Figure 4.10: Block-spin transfomiation is applied to a lattice repeatedly at T = 
0.99Tc. Each time the number of spins is reduced by a factor of 9, elucidating 
the behaviour of the system at a larger scale. The correlation length decreases 
under successive transformations. Clearly seen is the decrease of the correlation 
length and a suppression of the fluctuations. Under each transformation the system 
appears more ordered, the system flows under the renormalization transformations 
towards zero temperature. (After K.G. Wilson, Sci. Am. 241, 140 (1979)) 
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Figure 4.11: As Fig. 5.9 but for T = 1.22Tc. Under each transformation, the 
correlation length decreases, the system flows towards infinite temperature. (After 
K.G. Wilson, Sci. Am. 241, 140 (1979)) 
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Figure 4.12: As Fig. 5.9 but for T = T^: The correlation length is infinite 
and remains infinite under successive transformations, the system remains at the 
critical temperature. (After K.G. Wilson, Sci. Am. 241, 140 (1979)) 
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T=0 




f =00 



Figure 4.13: Block-spin renormalization at i/ = 0. Upon successive decimations 
of an initial configuration at T < Tc the system becomes more and more ordered 
as if the temperature were lower and lower. Eventually, the system becomes com- 
pletely ordered, this is the T = fixed point. Upon decimations of a state at 
T > Tc, the system becomes more and more disordered, the correlation length 
decreases as if the temperature were higher, until the correlation length becomes 
on the order of the lattice spacing. This is the infinite-temperature fixed point. 
A system at Tc stays, upon successive decimations, at Tc- This is another fixed 
point, the critical point. 

the properties in that region and then transform back to the original lattice. This 
is the basic idea of the renormalization-group methods and we shall exploit this 
idea on some simple examples. 

4.3.1 A Trivial Example: The d=l Ising Model 

In the renormalization-group methods the temperature changes under successive 
transformations, therefore it is convenient to work with the reduced Hamiltonian, 
we divide H by ksT. The reduced Hamiltonian of an Ising Hamiltonian with NN 
exchange interaction in an external field is 



n^-KY,SiSi^i-hY,Si, K^— h=—. (4.65) 
. . knT knT 



Recursion relations. If we want to make several successive renormalization- 
group (RG) transformations, the structure of Ti must not change after each trans- 
formation. We are thus seeking for such a transformation from the old Hamil- 
tonian H to the new Hamiltonian H' that keeps the structure of the Hamiltonian 
unchanged under successive transformations. In addition, the transformation must 
be such that it does not change the partition function. That is necessary because 
this ensures that the physical properties of the system are not changed. The parti- 
tion function of the one-dimensional Ising model in an external field is written in 



J 



H 
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— O— X— O— X— O— X— O— X— O— X— 
1 23456789 10 

Figure 4.14: Decimation on the d — 1 Ising model. 

the form 

Z{N,K,h) = 5^e-^ = 5:e^Si^^^^+i+'^S^^^ 

{S} {S} 
^ ^gJf(5i52+5253)+/i52+|(5i+53) 

{S} 

ygi^(S'3S4+54S5)+/l54+| (53+55) y . . . ^4 gg-J 

We divide the lattice into the sublattices x and o, denote the spins on the x sub- 
lattice by cTj (Fig. 4.14), 

2; ^ ^ ^gO-2[K(5i+53)+/i]+|(5i+53)ga4[Js:(53+55)+fe] + |(53+55) 
{S} {a} 

(4.67) 

and carry out a partial trace, i.e., we sum over 2^/^ configurations of the {a} spins 
sitting on the x sublattice: 

Z(N,K,h) = ^[e(^+l)(5i+53)+/. + e(-/f+t)(5i+53)-fe] x---. (4.68) 

{S} 

The entire partition function must not change and H must keep its structure, there- 
fore we have the condition 

Z{N,K,h) = e^^(-^''^)Z(-,ir',/i') =e^^^e-^' (4.69) 

2 {5} 

where 

n' = -K' ^ S,Si+2 -h' Si (4.70) 

ieodd ieodd 

and g is?L constant. This will ensure that the physical properties will not change 
under successive iterations. The two conditions (4.70) and (4.69) imply that 

^{K+\){Si+S,+2)+h _^ ^{-K+\){S,+S,+2)-h 
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^ ^K'SiSi+2+^(Si+Si+2)+2g _ (4.71) 

The last equation has to hold for any 5'^, Si^2 = ±1: 

^2K+2h ^-2K ^ ^K'+h'+2g 

^2K-2h _^ ^-2K ^ ^K'-h'+2g 

^+h^^-h ^ Q-K'+2g^ (4_72) 

from which we get: 
c — 



g2K+2h _|_ g-2ii" 
g2K-2/i _|„ Q-2K 



^8g ^ (^^2K+2h^^-2K^^(^^2K-2h^^-2K^^(^^h^^-hy (473) 



or: 



1, cosh(2ii' + /i) cosh(2/s: - /i) 
K — — In 



4 cosh h 

, 1, cosh(2X + /i) 

^ = ^+2^%osh(2X-/.) 

^ = ^ In f 16 cosh(2X + h) cosh(2X - h) cosh^ /il . (4.74) 

8 L -I 

Equations (4.73) or (4.74) are the recursion relations and determine the fixed 
points and the flow diagram of the system. 

In each iteration, the number of degrees of freedom is reduced by one half, the 
new Hamiltonian H has only one half of the previous spins, 

N' ^N/b 6 = 2, (4.75) 

and the lattice spacing is increased, 

a' = ba. (4.76) 

Other lengths which are measured in units of the lattice spacing are reduced by a 
factor b. For example, the new correlation length is 

r = (4.77) 

The remaining spins on the decimated lattice interact with their new nearest neigh- 
bours through the renormaUzed coupUng constants K' and are subject to renor- 
malized fields h'. 
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Fixed points. At a fixed point the parameters (the coupUng constants - in our 

case K and h) do not change under successive decimations. This means that the 
system at a fixed point stays at this point in the parameter space. With the new 
variables 

X = e-^^ y = e-^'^ z = e"^^ 

0<{x,y,z)<l (4.78) 

the recursion relations become 

+ 



X — X 



y' ^ y 



{x + y){l + xy) 
X + y 
1 + xy 



z' = z^xy^ ^— ^ (4.79) 

{x + y){l + xy){l + yy 

The recursion relations for x and y do not depend on z. Physically, this means 
that the singular behaviour of the free energy does not depend on a shift in the 
energy scale. Therefore, we first investigate the flow diagram and fixed points 
in the (x, y) plane, see Fig. 4.15. From the first equation it is easy to see that 
x' — X if x = 1, i.e., x^l is a fixed point of this system for arbitrary y . We 
denote the fixed points by a *. There is a line of fixed points at x* = 1 for 
< y* < I, corresponding to an infinite-temperature sink. Along this line, the 
spins are completely disordered and ^ — >^ 0. If we set x = 0, the second equation 
reduces to y' — y^, which has two fixed points, one at y* = 1 and one at y* — 0. 
The fixed point at {x*,y*) = (0, 1) (this means H — and T = 0) corresponds to 
the critical point where ^ ^ cxd. This is correct because we know that the d = 1 
Ising model has Tc = 0. This fixed point is unstable and any point in its vicinity 
will flow under RG transformations away from it towards the fixed points on the 
line at X* — 1. The critical behaviour of the d — 1 Ising model is governed by 
the least stable fixed point at {x*, y*) = (0, 1). The point at (x*, y*) = (0, 0) is 
another fixed point which corresponds to T = and H — oo. At this fixed point, 
the spins are completely ordered. 



4.3.2 RG Transformations on the d = 2 Ising Model 

The recursion relations of the d = 1 Ising model could be formulated exactly. 
Usually, however, farther-range interactions are introduced under successive dec- 
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Figure 4.15: Schematic flow diagram of the d = 1 Ising model. (After J.M. 
Yeomans, Statistical Mechanics of Phase Transitions, Oxford 1992.) 



imations and approximations have to be made in order to get the closed-form 
recursion relations. 

As an example, let us apply the decimation on the d = 2 Ising model on 
a square lattice in the absence of an external field. We start with NN (nearest 
neighbours) exchange interactions and sum the partition function over all possible 
configurations of the a spins (these are the spins that will be decimated, they are 
denoted by x in Fig. 4.16). 



{S} {S} {a} 

_ J~~J ^QK{Si+x+Si-a:+Si+y+Si-y) _j_ ^-K(Si+x+Si-x+Si+y+Si-y) 

{S} i 



(4.80) 



S±x and S±y are the "right/left" and "upper/lower" neighbours of Uj, respectively. 
The conditions are that under decimation Z must not change and that H should 
have the same structure as the original Hamiltonian. However, we soon realize 
that we cannot fulfill the above conditions with the renormalized form of the initial 
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Figure 4.16: Decimation on the d = 2 square lattice. The spins denoted by x 
are decimated. K is the NN interaction and L the NNN interaction. After deci- 
mation, the lattice constant increases by 6 = \/2, and previous NNN become NN 
(nearest neighbours). Solid lines represent the original lattice and dashed lines the 
decimated lattice. 

Hamiltonian. After decimation we get a more complicated Hamiltonian: 

T^' — — 5Z i-^' i^i+x^i+y + ^i+y^i-x + ^i-x^i-y + ^i-y^i+x) 

i 

~l~ L (^Si^x^i—x ~l~ ^i+y^i—y) ~l~ Q ^i+xSi+ySi—xSi—y + C ] , 

(4.81) 

where K' is the new NN interaction, L' the new NNN (next nearest neighbours) 
interaction, Q' the four-spin interaction, and C a constant. Only with the Hamil- 
tonian of this form, it is possible to satisfy the condition Z' — Z for arbitrary 
configurations of S±x,y. Although we started with only NN interactions, decima- 
tion generated NNN interaction and even four-spin interactions. So, we have to 
start from the beginning again with a more complicated Hamiltonian which in- 
cludes all these interactions. In the following we shall neglect Q' and assume that 
L is small. In this case we obtain the following approximate recursion relations: 

K' ^ L + 2K^ 

U ^ (4.82) 
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Figure 4.17: The renormalized interaction after a decimation is L + 2/^^. 



In fact, these relations could be guessed immediately, see Fig. 4. 17. The new near- 
est neighbours are the previous NNN, therefore the term L in the first equation. 
Besides, they were connected by two pairs of consecutive bonds, this brings 2K'^. 
In the second equation, the new next nearest neighbours (NNN) are connected by 
two consecutive NN bonds on the old lattice - each of them contributing a factor 
K. 

The recursion relations 4.82 have two trivial fixed points, {K*, L*) = (0, 0) 
and (K*,L*) = (oo, oo) and one non-trivial fixed point at (1/3, 1/9). The first 
one is the infinite-temperature, paramagnetic fixed point, the second one is the 
zero-temperature, ferromagnetic fixed point. The non-trivial fixed point corre- 
sponds to the critical point of the system. 



4.3.3 General Properties of the RG Transformations 

So far, we have seen that the renormalization-group transformations work by 
changing the length scale and by reducing the number of the degrees of freedom. 
Only at criticality the system remains unaltered under the change in scale and the 
critical behaviour is determined by a fixed point of the transformation. 

We start with an effective Hamiltonian, which is characterized by a set of 
coupling constants and fields K, L, ■ ■ ■. Each set of these constants defines a point 
in the "Hamiltonian space," the space of coupling constants. The renormalization- 
group transformations work in this space, 



w = nn. 



(4.83) 



The renormalization-group operator TZ reduces the number of degrees of freedom 
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from N to N'. This defines the scale factor b: 

N' 



(4.84) 



In the previous examples, the reduction was made in real space either by block- 
spin transformation or by decimation; later we shall do the reduction also by inte- 
grating out large momenta (short wavelengths) in momentum space. 

The essential condition to be satisfied by any RG transformation is that the 
partition function must not change. 



(4.85) 



Therefore the total free energy does not change but the free energy per unit cell 
(spin) increases as 



f{u')^b''f{n). 



(4.86) 



All lengths, measured in units of the new lattice spacing, are reduced by the factor 
b. Thus, the correlation length scales as 



(4.87) 



Fixed points. The iteration of IZ traces a trajectory in the Hamiltonian space, 
this is called the "Hamiltonian flow." The trajectories end in fixed points H*, de- 
fined by 

n* = nn* ^ n' = n = n*. (4.88) 

At the fixed points, the system is invariant under subsequent scale changes and 
thus the correlation length does not change 

= e = C- (4.89) 

Let us now take a closer look at the Hamiltonian flow near the fixed points. 
We introduce a reduced Hamiltonian which we write in the general form 

H = jl-$ (4.90) 

and investigate the flow in the parameter space spanned by the generalized fields 
jl, conjugate to the operator which contains product of operators (like spins). 



108 



I. Vilfan 



Statistical Mechanics 



Chapter 4 



• Example: In the case of the d = 2 Ising model on a square lattice, the 
components of cf) include the terms linear in S (the conjugate filed is h), 
quadratic in S (coming from interactions with NN, NNN neighbours), four- 
spin terms, etc. The corresponding conjugate fields are the magnetic field h 
and the coupling constants K, L, etc. In the beginning, many components 
of jl are zero, but under successive transformations, more and more complex 
interactions emerge and an increasing number of components of jl becomes 

Under successive RG transformations, the Hamiltonian evolves and the system 
moves through the parameter space: 

jl' = TZjl (4.91) 

until it reaches a fixed point when: 

ll'^ll=ir. (4.92) 

Every particular type of criticality corresponds to a particular fixed point and de- 
viations from the critical point are related to deviations from the corresponding 
fixed point. To get insight into the critical properties we must investigate the flow 
in the Hamiltonian space in the vicinity of the fixed points. 
Near a fixed point we linearize the operator TZ. We write 

jl = jT + 5jl 

jH = jT + Sjl' (4.93) 
where the small deviations from jl* are related by 

on 



5jl^jl* + jC(jr)5jl (4.94) 



6jl' = C{jr)Sjl. (4.95) 

The linear matrix C is evaluated at the fixed point and is thus a constant. The 
eigenvalues Aj and the eigenvectors i7i of C, 

= \Vi, (4.96) 

determine the critical properties of the Hamiltonian. 
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The eigenvalues Aj are functions of b. Two successive transformations with 
scale change b make a total scale change 6^, they are thus equivalent to a single 
transformation with the scale change 5^. 

Xi{b)Xi{b) = Ai(62) (4.97) 

and the eigenvalues must be of the form 

Xi(b) = (4.98) 

Xi are the critical exponents and, as we shall see, they are closely related to the 
usual critical exponents a, (3, ■ ■ ■. So, if we know Aj, we can obtain the critical 
exponents Xi from 

• Example 1: For the d = 1 Ising model, jii = x and 112 = y- After lineariza- 
tion of the recursion relations (4.73)around the fixed point at jl* = (0, 1) we 
get: 

Sx'^ASx Sy'^25y. (4.100) 

Hence, C is diagonal in this representation. The eigenvalues and corre- 
sponding eigenvectors are: 



Ai = 4 ui 



A2 = 2 i^i = ( ? ) • (4.101) 



1 

It follows that the two critical exponents are both positive: 

= 2 X2 = 1. (4.102) 

Example 2: For the = 2 Ising model, pi = K and 112 = L, and after 
linearization Eqs. (4.82) become 

5K' = AK*6K + 5L 

5L' = 2K*5K. (4.103) 
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For the non-trivial fixed point at (K*, L*) — (|, |), the RG transformation 
matrix is: 



(4.104) 



AK* 1 




4/3 1 


2K* 




2/3 



The eigenvalues of £ are: 

A, = . 1.72, 

3 

and the eigenvectors: 



1 



10 



2.58 
1 

-0.58 
1 



-0.387 (4.105) 



(4.106) 



(Question: what does negative A2 mean physically, how does the system 
behave under successive RG transformations if A is negative?) The corre- 
sponding critical exponents are: 



Xi 



^^ ln[(2 + v/lQ)/3] 



.In 1(2 



X2 



In 2 

10)/3| 



In 2 



57 



-2.74, 



(4.107) 



one is positive and the other one is negative. 



Flow in the parameter space. We first express Sjl near the fixed point at jT in 
terms of the eigenvectors Pi of C. In case of the d 1 Ising model this means that 



In general we write 6 jl as: 



5x' = Xi6x = iSx 
Sy' ^ X^Sy = 2Sy. 



(4.108) 
(4.109) 



The coefficients hi are the linear scaling fields. They are the components of the 
generalized field in the coordinate system of the eigenvectors. Under a renormal- 
ization 



(4.110) 
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where Aj = On the other hand: 



3ji' = 



(4.111) 



therefore 



(6>1). 



(4.112) 



It follows that the flow of the Hamiltonian in the parameter space close to the fixed 
point depends crucially on the exponents Xi. For positive Xi (Aj > 1), the scaling 
filed hi increases under repeated RG transformations and drives the system away 
from the fixed point. In this case, hi is a relevant scaling field. If the exponent 
Xi is negative, the corresponding scaling field hi decreases under subsequent RG 
transformations, thus driving the system towards the fixed point. Such fields are 
called irrelevant scaling fields. 

Thus the stability of the fixed point and the topology of the flow lines around 
it depend on the number of relevant and irrelevant scaling fields associated with 
it. To illustrate this. Fig. 4.18 shows the the trajectories in the parameter space 
with two irrelevant and one relevant scaling field. 

In general, any system which has initially the component of the relevant field 
different from zero, will after repeated RG transformations eventually be driven 
away from the fixed point (dotted lines in Fig. 4.18). Physically this corresponds 
to the system moving away from criticality. 

The points in which the relevant field vanishes, define the critical surface. 
All the points on the critical surface flow into the fixed point. The points on the 
critical surface represent systems with an infinite correlation length. Under RG 
transformations, the systems stays on the critical surface (^ remains infinite), but 
it moves towards the fixed point. 

• Example: The d = 1 Ising model has hi = x and /i2 = H- For the non-trivial 
fixed point at x* = and y* = I the critical exponents are xi = 2 > and 
a;2 = 1 > 0, that means that both, temperature and the external magnetic 
field, are relevant. 

In the case of the d = 2 Ising model on the square lattice, xi > and 
X2 < 0, as we have seen. It has one relevant scahng field (hi) and one irrel- 
evant field (h2). The relevant scaling field is related to reduced temperature 
whereas the irrelevant field is related to L, the coupling with NNN. Notice 
that here we considered the Ising model without external field. If the field 
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Figure 4.18: A fixed point with one relevant (h\) and two irrelevant (/i2 and /is) 
scaling fields. All the points on the two-dimensional critical surface flow to the 
fixed point (solid trajectories). The points slightly above (or below) the critical 
surface first flow towards the fixed point but eventually they are driven away from 
it as the relevant scaling field increases (dotted trajectories). 
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(c) 

Figure 4.19: Topology of the RG flow near fixed points with: (a) two relevant 
fields, (b) two irrelevant fields, and (c) one relevant and one irrelevant field. 

is included, then we have two relevant scaling fields and one irrelevant (the 
magnetic field is always relevant). 

Figure 4.19 shows possible scenarios in the case of two scaling fields. Ferro- 
magnets have two relevant fields, external field and temperature (or NN exchange 
interaction) and the fixed point is of the type shown in Fig. 4.19 (a). 

Universality. The critical behaviour is completely determined by the behaviour 
of the system close to the fixed point, where the transformation matrix can be 
linearized. The transformation matrix depends on the values of the scaling fields 
at the critical point and not on the initial values of the irrelevant scaling fields. 
Thus, all systems which flow close to the fixed point, independent of where in the 
parameter space they start, will exhibit the same critical exponents determined by 
the eigenvalues of the linear transformation matrix at that fixed point. The critical 
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exponents are universal and do not depend on the initial point in the parameter 
space (like the coupling strength). 



Critical exponent v. To illustrate how the critical exponent v is obtained from 
Xi, let us consider ferromagnetic systems. We know that there the critical point is 
specified by H and T. These are the relevant fields for magnetic systems. There- 
fore we identify hi by the reduced temperature, hi = t and /?.2 by the reduced 
magnetic field, h2 = h. Under RG transformations ^ and t scale as 

C'^i t'^b'^H (4.113) 

b 

The two equations tell us that ^^H is constant, it is invariant under RG transfor- 
mations. Therefore 

C(t)ocr^/^* (4.114) 
Since, by definition ofu,^(x t~^, we get 



1 

u = — 



(4.115) 



The last equation can be obtained also in a different way. Unless the system is 
initially at the critical point, the correlation length will decrease under successive 
block-spin RG transformations. The system will eventually reach a point where 
the correlation length ^' will become equal to the size of the renormalized unit cell, 
so that the nearest-neighbour spins become uncorrelated. This point corresponds 
to a certain temperature to which cannot depend on the initial temperature. Sup- 
pose, the system is initially at the temperature t, and arrives at after n iterations. 
Then, 

to^b^'^'H and ^' = 1 = =^ ^ = fe" (4.116) 

From these equations we get Eq. 4.1 15 again. 

• We return again to the d = 2 Ising model. Above, we considered the model 
at /i = 0. Thus, the exponent Xt corresponds to the relevant exponent xi = 
1.57 from Eq. (4.107) and therefore u — 1/xi — 0.64. This is better than 
the mean-field value of 1/2, but still far from the exact value v = \. The 
discrepancy has its origin in the approximations we have made to obtain the 
RG recursion relations. 
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Other critical exponents. The original and renormalized Hamiltonians Ti and 
Ti' correspond to positions in the parameter space, pi and jl' . We apply relation 
(4.86) to the singular part of the reduced free energy per spin: 

9s{il')^b''g,ifl). (4.117) 

Near a fixed point p, and p' can be expressed in terms of the scaling fields hi and 
h\ which are related by (4.112). Hence eq. (4.117) becomes: 

gs{liuh2,...) = b-^gs{b^'h,,b^'h2,...), (4.118) 

where gg is a homogeneous function of p. This is the scaling form of the singular 
part of the free energy. This equation holds for arbitrary b because we can repeat 
transformations many times. In fact, b can be any real number > 1 ! Again we 
identify hi — t, h2 — h and assume that other scaling fields are irrelevant. Then, 
eq. (4.118) becomes 

g,{t,h)^b-%{b^H,b^^h). (4.119) 
Now, to come to the specific-heat critical exponent a, we remind that 

CH = (^^) = gtt{h = 0) oc (4.120) 

Differentiating Eq. (4. 119) twice with respect to temperature and putting the scal- 
ing field h equal to zero gives: 

dr 

gu{t, /i = 0) « 6-Vr(r, 0)— = b'^+^^'grr, (4.121) 

where we have used r = b^^t. Initially, t is small whereas 6 > 1 so that the 
product b^^t increases with each iteration. We repeat the iterations until 

b'^'lt]^!. (4.122) 

In this way, the arguments of grr and therefore also the function g^r are constant, 
all the temperature dependence of the specific heat is in the prefactor b-^+'^^i ; 

gu{t, h = Q)(x |t|('^-2-i)/-i^,,(±l, 0) (4.123) 

and 

a = 2 - — . (4.124) 

Xi 
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Similarly we find 

7 = 



Xl 

S^-^ (4.125) 

a- X2 

We have expressed five critical exponents in terms of two independent expo- 
nents Xl and X2. That means that there must be relations between the exponents 
a - ■ - p. One can easily verify the scaling laws: 

a + 2/? + 7 = 2 
7 = /5(<5-l) 
2-a = diy. (4.126) 



4.4 Momentum-Space Renormalization 
Group 

RG methods in momentum space are much more powerful and more widely used 
than their counterparts in real space. The problem is only that they are much 
more complicated than the above real-space prototype examples. We shall define 
the RG transformation for the Landau model, that means that we shall now use 
continuous instead of discrete spin variables. The main reason is that discrete 
spins become essentially continuous after (block-spin) transformations. 
We start with the free energy of the form: 



7i(/i,m(r)) 



dV 



V 



g\Vm{r)\^ — m{r)h 



-|--rm^(r) -|- Mm'^(r) 



(4.127) 



This is also called the Landau-Ginzburg-Wilson (LGW) Hamiltonian. The main 
procedure of all the RG transformations is "coarse graining"(i.e., transformation 
to larger unit cells) followed by rescaUng (to make the system look like the original 
one). In real-space RG this was made either by introducing block spins or by 
decimation. Now we shall make coarse graining in momentum space. Instead 
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of expanding the size of the unit cell in real space, we shall shrink the size of 
the Brillouin zone in the momentum space. In the following the method of the 
momentum-space renormalization group will be demonstrated on the Gaussian 
model. 

4.4.1 The Gaussian Model 

The Hamiltonian of the Gaussian model in the momentum space is: 

T.i^ + 9q')\m{q)\'-hm{0). (4.128) 

|91<A 

We assume a homogeneous magnetic field and we are only interested in the stabil- 
ity of the m = solution (r > 0), therefore we take u — 0; this will simplify the 
summation over g). A is the cutoff momentum (Brillouin zone boundary). Taking 
u — means that the modes with different momenta become uncoupled. 

The corresponding partition function is equal to the functional integral over al 
functions m(r). In momentum space, the functional integral becomes: 

Dm(r)^ n / dm(q) (4.129) 

and the partition function: 

Y{t,h)= TT rdm(g)e-^^i^i<A^''+^^'^l"^('')l'-'^"^(°^ (4.130) 

The RG transformation consists of "integrating out" the wavevectors that are 
outside the sphere of the radius A/b and then rescaling. One RG cycle consists of 
three steps: 

1. Integration. After "integrating out" the momenta within a thin shell be- 
tween the radii A/b and A, the partition function (which should not change) 
becomes: 

y = e^ n /dm(g)e-^', (4.131) 

where the first term, cxp(Q), is a function of A, b, and coupling constants. It 
comes from the integration over momenta in the shell between A/b and A. 
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The second term is the partition function of the states with momenta \q\ < 
A/6. Notice that this separation into two contributions was only possible 
because we set u = Q, the partition function is then a product of independent 
Gaussian integrals. The new Hamiltonian, 

n'{m{q))^ E lir + gq')\m{q)\'-hm{0), (4.132) 

depends only on m{q) with |g| < A/6 and has the same structure as the 
original one. 

2. Rescaling. The second step restores the old cutoff by blowing up the radius 
of integration to the original value A in such a way that we change the 
variable of integration to 

q'^bq. (4.133) 

In this way the cutoff momentum is restored back to A. This step corre- 
sponds to putting a' — ha in real-space RG. The Hamiltonian now reads: 

n'{m)^h-<' E \{T + g'^^)\m{{)\'-hm{Q). (4.134) 

|9'|<A ^00 

3. Normalization. The step 3 restores the standard normalization of the order 
parameter. This step corresponds to putting the new spin equal to 1 after 
block-spin or decimation in real-space RG. Renormalization is restored with 
g' — gif we set: 

The transformed Hamiltonian and the partition function (which must not 
change) now look like 

n'{m)= E ^(r' + ^g'')|m'(g')r-/iW(0). (4.136) 

k'KA 

and 

y = e" n f ^rn{q')Q-^' . (4.137) 
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In this way we integrated out large momenta and transformed H to its orig- 
inal form. The recursion relations are: 

r' — b'^r 

Since 6 1 we can write 6 = 1 + dr, thus, 

b"" ^ 1 + nlnb = 1 + ndr. (4.139) 

We can say that r counts the number of RG cycles. Then, r and h obey the 
following differential equations: 

dr 
dr 
dh 
dr 

with the solution 



2r 

d + 2 



h 



(4.140) 



r = roe 
h = hoe'"'^ 



xt = 2 



d + 2 



(4.141) 



Since q is continuous (for macroscopic systems), the RG transformations was 
carried out in infinitesimal steps (6 k, 1) and the problem was reduced to solving 
a set of differential equations. 

The Gaussian model has two trivial fixed points at h* — 0, one at r* = oo 
(very high T, m = 0) and one at r* = —oo (very low T, where m ^ oo because 
we set M = 0!), and one nontrivial fixed point (the critical point) at h* = r* = 0. 
At the nontrivial fixed point the critical exponents are both positive and r and h are 
both relevant scaling fields. From Xi, we shall calculate other critical exponents 
in the Section 4.4.3. 

As we shall see in the next Section, for < 4 the Gaussian model is unsta- 
ble with respect to the parameter u coming from the four-spin interaction in the 
Hamiltonian. 
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4.4.2 The Landau- Wilson Model 

In the Gaussian model, we have neglected the fourth and higher-order terms in the 
Landau Hamiltonian. The next step is to include also the term. This will lead 
to non-MFA exponents, as we shall see. Here we shall bring only an outline of 
the RG procedure and discuss the results. We start with the LGW Hamiltonian in 
real space. 



/■ rl 1 
n{m) = / dV -c/|Vm(r)p + -rm^{r) + um^{r) 



(4.142) 



where u > ^. We first integrate out an infinitesrmally thin shell of the thickness 
(5gf in momentum space by taking 6 fa 1. Then, 

(5?/A = (5A/A = 51nA = 6-lRiln6 (4.143) 

and we write 

m(r) m(r) + 5m(r) (4.144) 

where 6m{r) contains the Fourier components to be integrated out and m the rest. 
After some approximations we get the following renormalization-group (differen- 
tial) equations for r and u [Huang, 1987]: 

^'^^^ = 2r + 12A'^-^.rl 



dr V A2. 

"'"^ = (4 - d)ii - SeA'^-V. (4.145) 



dr 

To analyse these equations, it is convenient to introduce the dimensionless 
coupling constants: 

r u 
^=A^' y=]^a (4.146) 



then, the RG equations become: 

dx 
dr 
dy 
dr 



2x + 12t/(l - x) 

ey- 367/2 (4.147) 
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where 



e = 4 - (/ 



(4.148) 



We shall assume that e is small. The system has two fixed points at finite x and y, 
fixed by dx/dr = dy/dr = 0: 



The Gaussian fixed point is at a; = |/ = whereas the "nontrivial" (also called 
the "Wilson-Fisher") fixed point lies in the upper {y > 0) half of the (x, y) plane 
only if e > 0, (i < 4. In the lower half of the {x, y) plane, u is negative and the 
system is unstable. The nontrivial fixed point approaches the Gaussian fixed point 
as e ^ 0. 

In the neighbourhood of the fixed points we linearize the RG equations 



X* = 
y* = 

X* = -e/6 
y* = e/36 



Gaussian fixed point 



'Nontrivial" fixed point. 



(4.149) 



X = X* + 5x, 



y = y* + 6y 



(4.150) 



and get the following linearized recursion relations: 



^dx = {2 - 12y*)5x + 12(1 - x*)dy 
^5y = ie-72y*)5y 



(4.151) 



which are written in the general form as: 




(4.152) 



At the Gaussian fixed point the linear matrix C is: 



2 12 
e 



(4.153) 



The matrix C has the following eigenvalues and eigenvectors: 
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X2 = e ^^°^(-(2-e)/12) ^^'^^^^ 

The corresponding scaling fields are hi and /i2- Under the (infinitesimal) RG 



transformations the scaling fields change as dhi/dr — Xihi 



(4.155) 



Thus, the critical exponents are equal to the eigenvalues Aj: = Ai = 2 and 

^2 = A2 = e 

In the neighbourhood of the nontrivial fixed point the linear transformation 
matrix C is: 



2(1 -6y*) 12(1 -x* 
e - 12y* 

It has the following eigenvalues and eigenvectors: 



2(1 -e/6) 12(1 + e/6) 
-e 



(4.156) 



Ai 
A2 



2 -e/3 



— e 



6+£ 



l+e/3 
1 



(4.157) 



Now we have the following picture. For d < 4 (e > 0) the nontrivial fixed 
point lies in the upper half of the {x, y) plane whereas the Gaussian fixed point 
is (as before) dX x — y — see Fig. 4.20. For the nontrivial fixed point, /i2 is 
irrelevant (it does not affect the critical exponents) and hi — r (x tthe relevant 
scaling field whereas for the Gaussian fixed point, both fields are relevant. Thus, 
the line z/2 is the critical line. For d> A{e < 0)the nontrivial fixed point lies in the 
lower half of the [x, y) plane, i.e., in the unphysical region {u < 0), see Fig. 4.21, 
and the critical behaviour is governed by the Gaussian fixed point, d — Ais the 
special case when both fixed points coincide. The linearized recursion relations 
give X2 — 0, /i2 is marginal, and X2 cannot tell us the direction of RG flow. In 
this case we must include higher-order terms which were neglected when we 
linearized the transformation matrix. These terms tell us that the scaling field /i2 
is irrelevant for u > at d = 4. 



4.4.3 Critical Exponents 

We have seen that the critical exponents of the nontrivial fixed point are 

xt = 2 -e/3 
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Figure 4.20: Fixed points and flow diagram of the Landau-Ginzburg-Wilson 
model for < 4. The Gaussian fixed point at (x*, y*) = (0, 0) has the mean- 
field exponents. The non-trivial fixed point at (x*,y*) = (— e/6, e/36) is unstable 
against the relevant perturbation 6x t and stable against the irrelevant perturba- 
tion 5y. The system at A has t < and it flows towards the low-temperature fixed 
point, r becomes increasingly negative. At the point C the system has t > and 
it flows to the high-temperature fixed point, r becomes increasingly positive. The 
point B is on the critical line and the system flows towards the non-trivial fixed 
point. (After K. Huang, Statistical Mechanics) 
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yocu 




Figure 4.21: Fixed points and flow diagram of the Landau-Ginzburg- Wilson 
model for d > A. The nontrivial fixed point is in the unphysical region. To 
stabilize the system in the lower half of the (x, y) plane, one needs an term. 
(After K. Huang, Statistical Mechanics) 



Xh = l + d/2. (4.158) 

In fact, we calculated Xh only for the Gaussian fixed point in the Eq. (4.138). 
However, the term with m'^ does not affect the value of Xh so that Xh of the non- 
trivial fixed point is the same. With these values of Xi we get from Eqs. (4.115, 
4.124, 4.125) to the first order in e for d < 4: 

e 
6 

1 _ e 

2 ~ 6 

3 + e 

1 e 

2 ^ 12 

0. (4.159) 



a = 

/9 = 

7 = 

6 = 

V = 

V = 
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For c? > 4 (e < 0) we only have the Gaussian fixed point in the physical region 
(u > 0) with the exponents 

xt ^ 2 

Xu — 6 

Xh = l + d/2. (4.160) 
With Eqs. (4.115, 4.124, 4.125) we find: 

. d 
" - 2-2 

7 = 1 



{S = 3 + e) 



1 

2 



r; = 0. (4.161) 

The above expressions for P and 5 are wrong, because we have neglected the 
term in considering the Gaussian fixed point. Below Tc, r is negative and for 
m to be finite, the term um^ with u > must be included in the Hamiltonian. 
Although u is an irrelevant field, it helps to stabilize mat a finite value and in this 
way influences the critical behaviour. By treating the um'^ term in the Gaussian 
approximation (independent modes) we get the usual mean-field values for the 
exponents /3 and 77: 

5^3 (4.162) 

for all d > A. These results are correct and also agree with the results of the 
Landau (MFA) model in Sec. 4.1.4. 

The Gaussian model gives the same critical exponents as the mean-field theory 
except for the specific-heat exponent a. Therefore the Gaussian fixed point is 
identified as the mean-field fixed point. In the Gaussian model the specific heat 
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diverges (a > 0) for d < 4. This divergence is caused by the long-wavelength 
fluctuations which were neglected in the calculation of the free energy in the MF 
approximation. 

In d = 3 (e = 1) we find a 0.17 and u 0.58. The results of this paragraph 
are based on an expansion of d around d — A and are valid for small e > or for 
d > A. Nevertheless, we put e = 1 when d = 3. That means that e is not small and 
higher-order terms in e must be included. This is then the e— expansion, which 
gives good values for the critical exponents. 
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